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Abstract 

Following the discovery that asteroid (3753) Cruithne was a coorbital companion of the Earth, a new theory of 
coorbital motion has been developed whereby planets or satellites can maintain companion objects in the same 
orbit as themselves. This has led to the prediction of hitherto unknown types of stable motion, all of which are 
seen in the evolution of specific near-Earth asteroids. The slow diffusion of such objects through the Earth's 
coorbital region is shown to lead to temporary capture, suggesting the existence of undiscovered retrograde 
Q^ ' moons of the Earth. 
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1. Introduction 



Lagrange^ demonstrated the existence of five equilibrium positions in the three-body problem. Those that form 
an equilateral triangle with the Sun and the planet, the leading L4 point and trailing L5 point, are stable to 
small displacements. However, it was not until the subsequent discovery of the first Trojan asteroid, (588) 
Achilles, near the L4 point in the Sun- Jupiter system that a real example of such motion was first observed. 
There are now known to be at least 400 Trojan asteroids in the orbit of Jupiter, and at least one, (5261) Eureka, 
\^ • near Mars's L5 point;^ there are also three examples of Trojan moons in the Saturnian system.^ In fact, all the 
known Trojan objects are librating about L4 or L5 points and their orbits are often referred to as tadpoles (or 
' T orbits) because of their shape with respect to the equilibrium position, as viewed in the frame rotating with 
\ the mean angular velocity (mean motion) of the planet (see Fig 1). 

When the angular amplitudes of the leading and trailing tadpoles are sufficiently large, the orbits merge near 
the unstable L3 point located 180° from the planet; the resulting paths in the rotating frame are referred to as 
horseshoes (or H orbits). A variation on such a structure is seen in the orbits of Janus and Epimetheus, the 
coorbital satellites of Saturn.'' In this particular case, a good theoretical understanding of the peculiar orbital 
dynamics is possible because of the low eccentricities (e < 0.01) and inclinations (/ < 0.4°) of both objects. In 
O \ contrast, it has recently been discovered that asteroid (3753) Cruithne, previously designated 1986TO, performs 
• a temporary horseshoe-like orbit with respect to the Earth. ^ This is a more extreme example that has been 
^ \ difficult to incorporate in an analytical theory because of Cruithne's large eccentricity (e — 0.515) and inclination 
(/=19.8°). 

In this article, we present aspects of the general theory of coorbital motion that accounts for all types of 
■ coorbital behavior.^ This shows that asteroid (3753) Cruithne is only a member of a more general class of objects 
H \ that are captured in the coorbital regions of the planets, and that involve new types of coorbital motion at 
large eccentricity and inclination. We also report the identification of these new orbits in the motion of specific 
asteroids. For the cases we have examined, the additional planetary perturbations can cause near-Earth objects 
to be trapped into the coorbital regions of the terrestrial planets without any risk of collisions. The structures 
we found provide a new paradigm for coorbital motion in the three-body problem and we believe that they have 
played an important role in such diverse problems as the origin of retrograde outer satellites and the accretional 
processes in the early solar system. 



2. Types of motion 

In the Sun-planet-asteroid three-body problem, the motion of an asteroid on a small eccentricity and inclination 
orbit is constrained by the value of the Jacobi constant through the 'excluded' regions associated with the zero 
velocity curves'" (Fig 1). The T and H librations are found to be limited to the annulus between the Lagrangian 
points Li and L2 defined as the coorbital region of the planet; its half-width equals the radius of the planet's 
sphere of influence: 

Ap = [mp/3(mp + ms)]^^^ ap (1) 
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where ap is the semi-major axis of the planet's orbit, and mp and mg are the masses of the planet and the 
Sun respectively. The outer boundary of the coorbital region is made up of irregular orbits and separates the 
H orbits from the passing orbits (or P orbits). 

For large eccentricities, the previous picture of coorbital motion is modified substantially because the Jacobi 
constant can no longer exclude spatial regions for the motion of the asteroid. In addition, the large eccentricity 
of the asteroid's orbit can lead to physical collisions with the planet. There are two facts that allows us to 
develop a theory that is valid for all types of coorbital motion: the first is the realisation that the asteroid's 
motion can be described as the superposition of two types of motion: (i) a fast three dimensional gyration with 
amplitudes of O(aAe) and 0(aAsin7) and a frequency equal to the asteroid's mean motion and (ii) the slow 
evolution of the mean position of the asteroid with respect to the planet. The latter is referred to as the motion 
of the guiding center and is represented by the conjugate variables: relative semi-major axis, Or = (oa — ap)/ap, 
and relative mean longitude, Ar = Aa — Ap . 

The second fact is the regularity of e, J and the argument of perihelion uj away from the boundary of the 
coorbital region; this allows us to derive the motion of the guiding center by perturbation analysis. The Jacobi 
constant can then be used to constrain the motion of the guiding center instead of the full motion, by averaging 
with respect to the fast gyration represented by Aa. The guiding center ar(Ar) can be written as: 

^2 0- S{K,e,I,c.) (2) 

3(mp -I- msj 

where C is a constant that defines the nature of the orbit, and S is related to the averaged interaction potential 
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S=—J (|rA - rp|"^ - FA • rp) dAA 



(3) 

ar=0 

In equation (3), ta and rp arc the position vectors, and the second term is due to the motion of the Sun; the 
distances are scaled to ap. The motion along the guiding center can now be viewed as that of a particle in the 
unidimensional potential well 5 (Fig 2). 

The fact that the shape of the guiding center depends on e, / and u, has three major consequences for 
the coorbital motion: firstly, the effective Lagrangian points L4 and L5 can be displaced appreciably from the 
equilateral configuration (Fig 2 b,c,d). In addition, they no longer correspond necessarily to the same Jacobi 
constant, which fact implies an asymmetry in the distribution of T orbits.^ 

Secondly, retrograde satellite orbits (or RS orbits) appear outside the planet's sphere of infiuence as they are 
the only family of bounded orbits around the planet that remain stable with increasing eccentricity^ (Compare 
a and b in Fig 2) . For the planar motion, it can be seen that this family is separated from the H and T orbits by 
the collision singularity, fa = rp, at Ar c±. 2e for e < 0.7. This implies that RS orbits and H orbits merge inside 
the inner boundary of the P orbit domain and therefore, that coUisional orbits are unstable in two dimensions. 

In three dimensions, however, physical collisions occur only when the nodes of the asteroid's orbit cross the 
orbit of the planet; this corresponds to a specific argument of perihelion given by cosw ~ e. Consequently, the 
collision singularity of S is generally removed and replaced by maxima whose existence and relative magnitudes 
depend on e, / and uj. 

This leads to the third major consequence of large e and /: the appearance of compound and transition 
orbits (Fig 2 c,d). The first correspond to the merger of H or T orbits with RS orbits; asteroid (3753) Cruithne^ 
is the first example of a compound H-RS orbit in the solar system. The second orbits correspond to the maxima 
of S; however, they are not long-lived because of the secular variations of e, / and uj, in time which modify the 
magnitudes of the corresponding maxima. Transition orbits therefore appear in the orbital evolution only to 
permit the passage between different orbit families. 

The secular stability that was found in the three-body problem,® together with the presence of asteroid 
Cruithne on an H-RS orbit, suggested the robustness of this new dynamics when exposed to planetary pertur- 
bations. This motivated our search for further examples in the solar system. 



3. New coorbital asteroids 

We have conducted a search in the catalogue of near-Earth objects (NEOs) which is maintained in the Minor 
Planet Center web-site^° for objects with semi-major axes closer than approximately Ap from the semi-major 
axis of a terrestrial planet. We carefully avoided objects with poorly determined orbital elements, choosing to 
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include in our sample only multi-apparition asteroids. We stress the need here for further observations that 
refine the orbits of the known near-Earth asteroids. Our selection criteria yielded a total of five objects, ~ 1% of 
the known NEO population, as favorite candidates for transient coorbital behavior. These asteroids are (3753) 
Cruithnc, (3362) Khufu, 1989 VA, 1993 WD and 1994 TF2. Hero, we present results for three of those objects 
(see Table 1) as representative of what is observed for the full sample. 

The orbits were integrated numerically^^ with the eight major planets (Mercury to Neptune) using a Runge- 
Kutta-Nystrom 12th order scheme with adaptive step-size control. This presents an advantage over the scheme 
used in ref 5 due to its ability to model close encounters, a vital part of the mechanism which gives rise to stable 
orbit transitions. The integration time-span was set to be 200,000 years centered on the present as a reasonable 
trade-off between integrator accuracy, which may deteriorate for longer time-spans, and the richness of the 
coorbital dynamics that we expect to observe. 

All the objects examined exhibited transient coorbital behavior. Figures 3 and 4 show the orbits' evolution 
during the coorbital capture of (3753) Cruithne and (3362) Khufu by the Earth, and 1989 VA by Venus. We 
confirm the existence of the new structures discussed in the previous section, such as compound orbits (Cruithne, 
1989 VA), displaced tadpole librations (Cruithne), retrograde satellites of the Earth (Khufu, Cruithne) and the 
transitions between different coorbital modes. Owing to the intervening periods of stochastic drift, mainly in 
Or, we cannot claim that we observe the actual evolution of the asteroids in the integration time-span. What 
we can say, however, is that each one of those asteroids is likely to have been or will be captured in such orbits. 
This claim is supported by our numerical experiments with slightly modified orbits for the asteroids (i.e. clone 
asteroids) . 



4. Discussion 

Unlike (3753) Cruithne which is currently on an H-RS orbit, well inside the coorbital region of the Earth, 
the remaining objects that we examined are not currently in coorbital libration with any planet. This fact, 
however, did not induce distinct types of global evolution in the sample. For instance, during the evolution, 
Cruithne happened to exit the coorbital region; in contrast, the other objects could be captured by the Earth 
or Venus. This shows that Cruithne's status is less special than previously thought: it is simply a member of 
the population of NEOs that evolve inside the inner solar system with a slowly varying semi-major ax;is."^^ 

This seemingly disappointing conclusion about the true nature of Cruithne's evolution can be viewed pos- 
itively in the context of the protection of the Earth from collisions with nearby asteroids. The evolution of 
the mentioned asteroids and their clones indeed suggests that the Earth and Venus are well protected from the 
class of NEOs that drift (in semi-major axis) towards and get close to the coorbital regions of the terrestrial 
planets. A key factor in this protection is the quasi-absence of the chaotic layer at the edges of the coorbital 
region, found at small e and /, and responsible for the instability during the crossing of coorbital region. ^'^^ 
This suggests the presence of undiscovered coorbital asteroids of the terrestrial planets, that display one of the 
new coorbital modes. In particular, the connection between retrograde satellite orbits and the classic types 
of motion is remarkable because it allows a fraction of the drifting NEOs to become, temporarily, retrograde 
satellites of the terrestrial planets. Clear examples of this behavior are the past trapping of Khufu that lasted 
35,000 years and the future 3,000 year RS phase of Cruithne's evolution. This same connection also offers a 
new route for the capture of the known retrograde satellites of the Jovian planets: the dissipative processes of 
early planetary formation such as the planet's growth^^ or the drag due to the circumplanetary nebula,^^ can 
produce permanent capture by reducing the Jacobi constant of a Trojan asteroid that undergoes a transient 
RS motion. This results in a distant retrograde satellite with large planetocentric eccentricity and inclination, 
which is reminiscent of the orbits of the known retrograde satellites of the Jovian planets.^'' 

The evolution of our sample shows that the coorbital asteroids of Earth and Venus with large e and / are 
generally expelled from the coorbital regions because of the close encoimtcrs with the other planets. In that 
respect. Mercury is a better candidate for finding permanent coorbitals because of its distant location from 
Venus. The presence of a specific population of minor bodies close to Mercury has already been proposed in the 
context of the vulcanoid hypothesis for the chronology of the geological evolution the planet.^* However, the 
dynamical constraints of the model assumed that the vulcanoids had small eccentricities and inclinations and 
the recent observational searches^^ produced negative results. In view of (1) the large eccentricity of Mercury 
(e = 0.2), (2) the action of nearby secular resonances and (3) the stability of the structures we found, the 
putative coorbital vulcanoids are likely to have stable large eccentricity and inclination orbits. The evidence of 
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such objects would argue for the early presence and slow depletion of a larger reservoir of bodies that contributed 
significantly to the cratering history of Mercury. 
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Figure captions 



1. Representative zero velocity curves for three values of the Jacobi constant and the mass ratio top/(top + 
ms) = 0.01. The locations of the Lagrangian equilibrium points L1-L5 are indicated by small open circles. 
The dashed line denotes a circle of radius equal to the planet's semi-major axis. The letters T (tadpole), 
H (horseshoe) and P (passing) denote the type of orbit associated with the curves. The regions enclosed 
by each curve (shaded) are excluded from the spatial motion of the asteroid that has the corresponding 
value of the Jacobi constant. Note that the curve between the Li and L2 points permits satellite orbits 
around the planet inside the Hill sphere of radius Ap (Eq 1). 

2. Types of coorbital motion. 5* is viewed as an effective potential well and the possible guiding centers are 
given by the levels 3(mp + m,s)C /Sm,p (light linos). The dashed levels correspond to transition orbits. 
T, H, RS and P denote respectively tadpole, horseshoe and retrograde satellite and passing orbits. The 
hyphenated designations denote compound orbits. Li and L2 are located at Ar — 0° and L3 at ±180°. The 
bold ticks denote the locations of L4 and L5 . The comparison of a (circular 2D orbits) and b (eccentric 
2D orbits, e = 0.3) shows the appearance of retrograde satellite orbits outside the Hill sphere and the 
displacement of L4 and is form the equilateral configuration Ar = ±60°. In 3D, (c: e = 0.3, / = 20°, 
u = 60°; d: e = 0.5, / = 30°, lu = 0°), orbit transition and compound orbits become possible. Note the 
asymmetry in the location and the maxima of L4 and L5. Because S{Xr,e, I, —oj) = S{—Xr,e, I,uj), the 
coorbital motion also includes orbits that are symmetric to those in c with respect to Ap = 0. 

3. Coorbital capture of near-Earth asteroids. The evolution of the relative semi-major axis, a^, and the 

difference in mean longitudes between the object and the relevant planet, Aj., are shown as functions of 
time. The plots correspond to (3753) Cruithne (a, b) and (3362) Khufu (c, d) coorbiting with the Earth, 
and 1989 VA (e, f) coorbiting with Venus. The asteroids are seen to transit between many types of 
coorbital motion (Fig 2) for time-scales ranging from thousands to tens of thousands of years. Note that 
in all three cases, the transitions occur without the large jumps in semi-major axis which are indicative 
of unstable close approaches. 

4. Coorbital librations of the guiding center. The trajectories observed in Fig 4 are shown in the {Xr,ar) 
plane, a: (3753) Cruithne, b: (3362) Khufu and c: 1989 VA. Note that the proximity of Khufu's guiding 
center (~ 1°) to the Earth (of radius ~ 9") does not result in a near collision because the full motion 
of the asteroid includes the fast gyration that was averaged out and whose smallest amplitude is ~ a as. 
This means that the asteroid is 0.5 AU away from the Earth. The temporary proximity to A^ = 0° 
implies that the trajectory associated with the fast gyration is centered on the Earth and does not precess 
secularly. 



Table 1: Osculating elements for three coorbital candidates at JD 2451000.5 as given by the JPL HORIZONS 
ephemeris generation system. oa, e, /, uj and 17 denote respectively the semi-major axis, the eccentricity, the 
inclination, the argument of perihelion and the longitude of the ascending node, a^ap is the difference between 
the semi-major axes of the coorbital planet and the asteroid; it is given in terms of the width of the planet's 
coorbital region Ap. OCC^^ (Orbit Condition Code) provides a rating of the quality of the computed orbit 
based on the available observations. In this scale is the best and 7 the worst with 1 the typical value for 
numbered asteroids. 
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